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Experience an electric force -eE directed away from the surface, ail escape, helping 
to discharge it. in Fig. 1(b), the spacecraft surface is nearly parallel to 'S, and 
almost all of the emitted electrons return to it., even though they still experience 
ah electric force directed away from it, -.These electrons therefore are unable to 
help discharge it, so a surface nearly parallel to t i s more likely to charge to a 
large negative voltage. Note that the component of E which is perpendicular «-o t 
results only in an-J-xfi drift parallel- to the surface. 

For any object much larger than 13 cm, the escape of secondary electrons will be 
strongly affected by this process. For example, most surfaces on the Shuttle are 
effectively "infinite planes" by this criterion. On the other hand, the average 
gyroradius of high-energy auroral electrons is comparable to Shuttle dimensions, so 
the deposition of these electrons onto Shuttle surfaces is likely to be only moder- 
ately inhibited. 

For a larger object (size >> 8 m) , deposition- of auroral electrons will also 
become strongly orientation-dependent, with both collection and escape of electrons 
now being inhibited on surfaces nearly parallel to B, This suggests that 
high-voltage charging of such surfaces may be more likely on objects of intermediate 
size than on either larger or smaller ones. In th- calculation Of Parks and Katz^, 
Katz and Parks , the tendency toward high-voltage charging increased with spacecraft 
size because in their model, ion collection increased less rapidly with spacecraft 
size than did electron collection. To determine which of these two effects predomin- 
ates will require more detailed calculations than have been done so far. 


As already mentioned, strong ion flow effects also are generally present in low 
orbit; the ion speed ratios (flow speed/most probable ion thermal speed) for H + at 
1 keV, H at 0.1 eV, and 0+ at 0.1 eV are 0.02, 1-8, and 7.3, respectively. Whenever 
the latter is the predominant ion species, ion collection on downstream surfaces will 
therefore be^srrongly inhibited. If a surface is simultaneously downstream and nearly 
parallel to B, as is likely to be the case in the auroral zones, then the tendency 
for high-voltage charging to occur on it will be greatly increased (Fig. 2) . 

To "straightforwardly" include B effects on secondary electron emission in a 
large two or three dimensional simulation program would involve the numerical integr- 
ation of very large numbers of secondary-electron orbits. The resulting computing costs 
usually would be formidable , especially since these orbits would have relatively large 
curvatures. A desirable alternative is to "parameterize" the situation by treating 
in advance a simplified but still sufficiently realistic model problem. In order to 
do this, we make the approximations described in the next Section. 

2. THEORY FOR £ NORMAL TO SURFACE 

We assume t^at th^ spacecraft surface is an infinite plane, and the electric and 
magnetic fields E and B outside it are uniform. In the work presented here, we also 
assume that the electric force -eE on electrons is directed along the outward normal 
to the surface; here e is the magnitude of the elementary charge. This assumption is 
to be relaxed in a later paper (J.G. Laframboise, to be published) in order to permit 
variations of potential along the surface to be taken into account. We assume that 
the secondary electrons are emitted with a Maxwellian distribution corresponding to a 
temperature T. The ratio i = l/l Q of escaping to emitted flux is then a function of 
two parameters: the angle 9 between the surface normal and the direction of B (Fig. 3), 
and a parameter describing the strength of E. A convenient choice for this parameter 
is the difference in potential across a mean secondary-electron gyroradius a = ( 1/eB) 
(irmkT/2) , divided by kT/e, where m is electron mass and k is Boltzmann's constant. 
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This quotient is: 


( 2 . 1 ) 


E f urn 

e - B \| 2kT 

where E ~ |ti| and B = |^| . 

This quantity also has an alternative, more useful interpretation: it is the- 

ratio of the magnitude |hxti|/B 2 of the Ex B drift speed, to one "ha If the mean 
thermal speed (GRT/irm)* 4 of the emitted electrons. It is useful to estimate the value 
of e for a high-voltage spacecraft sheath in low-orbit conditions. To do this, we 
use the sheath solution of Al'pert et al (Ref. 3, Table XXIV and Fig. 72). For a 1 kV 

and a 5 kV sheath around a sphere of radius 3m in a collisionless plasma having an 

ambient ion temperature of 0.1V, number density of 3 x 10 J cm” , and resultant (ion) 
Debye length of 0.43 cm, their results give*, respectively, sheath thicknesses of 2.6 
and 6.1 m, and surface electric fields E ~ 0.86 and 2.9 kV/m. Using B * 4.4 x 10” T and 

T= 3 eV for secondary electrons, we then obtain e - 33 . 9 and 114.2. Both of these 

are relatively large values, whose significance can be understood if we consider what 
would happen if e were infinite. 

In this limit, it is easy to show that secondary electrons would all escape 
unless B were exactly parallel to the surface (0 were 90°). This can be shown as 
follows. In this limit, secondary electrons would have no "thermal" motion. The 
(y,z) projection of their motion would then be similar to that shown in Fig. 4. This 
motion would be the sum of: (i) an-E-xB drift in the y direction (ii) a uniform 

acceleration along B, whose projection in the (y,z) plane would be upward (iii) ]ust 
enough gyromotion to produce a cycloidal path when combined with (i) , so that in the 
absence of (ii) , the electron would (just) return to the surface at the end of each 
gyroperiod. In the presence of (ii) , these "return points" are displaced upward by 
progressively increasing amounts (Fig. 4) , so the electron can never return to the 
surface, unless t is exactly parallel to the surface, so that the-upward component of 
-el: along $ vanishes. 

This result suggests that for large finite values of e (including the values 
calculated above) , electron escape is likely to be almost complete except for 0 very 
near 90°, where it should drop to zero very steeply. The occurrence of high-voltage 
charging in marginal circumstances may therefore depend very strongly on the precise 
orientation of a surface. 

The escaping secondary-electron flux is given by: 

1 = /J/f <v 0 )H(v 0 )v 0Z d3 ^o 

where: v 0 is the initial velocity of an emitted electron, f (v Q ) = d^n/d^v c is the 

velocity distribution of emitted electrons, n is a reference number density, and 

H(v 0 ) is equal to 1 for 'scaping electrons and 0 for those which retu rn to . the — 

surface. The emitted flux is: 

I * n (kT/2itm) **. (2 

o 

We also introduce the dimensionless Velocity: 

u » v(m/2kT) li . ^ 
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( 2 . 2 ) 
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Equation (2*2) then becomes: 
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which is in a form suitable for numerical summation. The quantities Uj. i # u i in 2' 
. . . . ,u. . are the values of u oz for which H changes between 0 and I ' for eAch 

u ox ana Uoy1 ,ax These values must be found by numerically determining which particle 
orbits reimpact the surface. These orbits can, however, be determined in analytic 
form, with time as a parameter. To do this, we use the coordinate system shown in 
Fig. 3, together with a y-axis (not shown) directed into the plane of the Figure. 
The equation of motion for an electron is: 



a +•*■ •> 

- (E + v x B) . 
m 


( 2 . 6 ) 


We solve this with the initial conditions £ = y=n =0, v^-v Q ^, Vy = v Q y, and v^ - v 0rj . 
We introduce the dimensionless variables : — J 


e 


x 





etc; 


x =x/a, $ = y/a, etc; 
T = = (eB/m) t . 


(2.7) 


In the present work, e x and Gy are both zero, but for future use. We have 

retained these quantities in the formulas below. We obtain: 


u _ = u sin 0 + u cos 0 ; 
oE, ox oz 

u on “ -Uox cos 9 + u oz sin 0 ; 
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( 2 . 8 ) 


Equations (2.8) can also be differentiated to find dz/dt . The numerical procedure 
for finding the quantities in Eq. (2.5) then involves calculating z and dz/dt 

at a succession of points along # an orbit (the electron will reimpact during the first 
gyroperiod 0 <t < 2tt if at all, So this interval always suffices), and making the 
appropriate tests on these quantities to find out whether the orbit reimpacts or 
escapes. For each u Q . and u Q y # j, this is done for a succession of values of u oz . 
These tests also yielci'ihe local minimum of £(t) if one exists. Whenever a change 
occurs between nc escape and escape from one such value of u oz to the next, an 
interpolation using these minima can be used to provide the corresponding value of 
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u Jim );* In ca5es where they are unavailable, the arithmetic mean of the two succes- 
sive 'u Q?< values is used. This completes the definition of the procedure used for cal-" 
culating the ratio I/I c of escaping to emitted flux. 


3. RESULTS AND DISCUSSION 

Escaping secondary-electron current densities, computed as described in Sec. 2, 
are shown in Table 1 and Fig,. 5. Each value of i - T/I Q was calculated using 191000 
orbits, evenly spaced in the intervals -4, 5 £ u ox i 4. Si, -4.5 1 u oy S 4.5, and 0 £ 

«on - 4.5, with points on the orbits calculated at intervals At~it/ 45. For 8 values 
each of e, and 0, the resulting calculation took 03 hr total on a Hewlett-Packard 
1000F minicomputer with Vector Instruction Set. The results are accurate to within 
about 0.5% or better. The result for e - 0 is just the analytic result i - cos 0 . To 
see why this is so, wc consider the electron orbit shown in Fig. 6, which has been 
fictitiously extended so as to pass through the surface and re-emerge from it. In 
the absence of an electric field (e * 0) , this orbit has the same speed at the 
re-emergence point C as at the- emission point A. Since we have also assumed that the 
emitted velocity distribution is isotropic, and therefore a function of speed only, 
the real orbit, for which C is the emission point, must carry the same population as 
would- the fictitious re-emerged orbit. The flux crossing the reference surface DE, 
which is±B, is therefore the same as if such passagesand re-emergences actually 
occurred, and is the same as if another reference surface FG, also J_B, were emitting 
electrons having the same velocity distribution. However, in reality, the electrons 
come from the real surface HJ, which is not J_ B, and all the electron-orbit guiding 
centers which are inside any given magnetic-flux tube through DE will also be inside 
the projection of the same flux tube onto HJ, and the ratio of the intersection areas 
of this tube with HJ and DE is just sec 0 . The ratio of escaping to emitted flux 
must therefore be the reciprocal of this, or cos 0 , as stated above. 

Also evident in Fig. 5 is the fact, mentioned in Sec. 2, that when e is large 
enough, electron escape becomes essentially complete except when 0 is very nearly 90°. 
In a real situation, E would not be uniform, but would decrease with distance from 
the surface, contrary to our assumptions. Our results can therefore be expected to 
overestimate electron escape. This would probably not be a large effect, but this 
presumption remains to be verified. An approximate compensation for it can be made 
by calculating £ using an electric field value which is averaged over the first mean 
gyroradius distance from the surface. 


The results in Table 1 are approximated to within 2.5% of I 0 by the empirical 
formula: 

a = 1 + l^t 1 * 1394 exp^O. 083725 -|l + tanh [1.9732 an ( 3^ 

- 0.07825 an [l + (e/8.5) 1 • 781 40 ]j ; 

1 + tanh [l. 49 an ( 

c » an (9O°/0); (3.1) 

i = cos [90° exp(-ac - be 2 )] . 

This formula also has the correct limiting behavior when e+0 or ®, or 0^-0° or 90° 
An approximation formula for the emitted flux is also available [Eqs. (5) and 



b = 0.38033e 0,85892 


exp ^2. 0988 ^ 
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(6) of Laframboisq ct al, Ref. 4, and Laframboise and Kami t sunk*, Rof. 5], 


4. 


CALCULATION OF SECONDARY “ELECTRON DENSITIES 
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Values of the ratio =I/I 0 of escaping to emitted flux, tar various values of 0 
he angle (xn degrees) between the surface norma) am\ magnetic field direction , 
and e, the nondime. isional repelling electric t strength. These two quantities 
appear in the table as THETA and EPS, respectively. These results are accurate to 
within about 0.5% or better; thus the differences between .999 and 1.000 in the Table 
are not significant. For 6 = 0°, i = l for all values of e. 
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figure 1. Effect of surface orientation on escape of emitted electrons. In (a), the 
spacecraft surface is perpendicular to the magnetic field B, and the emitted electron* 

.!* PerienCe an electric fo *' ce directed away from the surface, all escape. In 
(b), the spacecraft surface is nearly parallel to B, and almost all of the emitted 
electrons return to the surface, even though they still experience an electric force 
lrected away from it. Note that the component of E perpendicular to I results only 
m an ExB dritt parallel to the surface. 
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Figure 2. |Pacecraft simultaneously in a collisionless ion flow and a magnetic field 
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Figure 3. coordinate system for calculating electron escape £3u*e~ 
(not shown) is directed into the plane of the Figure. 


The y-coordinate 









